Schrodinger-Newton equation as a possible generator of quantum state reduction 
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It has been suggested by Diosi and Penrose that the occurrence of quantum state reduction 
in macroscopic objects is related to a manifestation of gravitational effects in quantum mechanics. 
Although within Penrose's framework the dynamics of the quantum state reduction is not prescribed, 
it was suggested that the so called Schrodinger-Newton equation can be used to at least identify the 
resulting classical end states. Here we analyze the extent to which the Schrodinger-Newton equation 
can be used as a model to generate a full, time dependent description of the quantum state reduction 
process. We find that when supplied with an imaginary gravitational potential, the Schrodinger- 
Newton equation offers a rationalisation for some of the hitherto unexplained characteristics of 
quantum state reduction. The description remains incomplete however, because it is unclear how 
to fully recover Born's rule. 



I: INTRODUCTION 

More than a decade ago Penrose and, independently, 
Diosi noticed that the inherent incompatibility of quan- 
tum mechanics' unitary time evolution with general rel- 
ativity's diffeomorphism invariance could have measur- 
able consequences on mesoscopic scales yj, [2| . Although 
the precise way in which general relativity would alter 
quantum mechanics on these scales is not known theo- 
retically, and proves to be extremely hard to access ex- 
perimentally 0, 0, @ , Penrose and Diosi use very general 
arguments to suggest that the process might be used to 
clarify the mystery of quantum state reduction. In a 
nutshell, the main argument is that gravity necessarily 
induces some uncertainty in the total energy of a quan- 
tum superposition state and thus allows quantum super- 
positions only a finite lifetime given by the inverse of 
that energy uncertainty. In turn, this can in principle 
lead to testable predictions about the timescale at which 
quantum state reduction is expected to occur [l], Q ■ On 
the other hand, this argument does not provide informa- 
tion on the dynamics of the reduction process itself. The 
equations that describe the transition from a well defined 
microscopic superposition state to a collapsed macro- 
scopic measuring apparatus remain unknown. Progress 
in this direction was made by proposing to use the so 
called Schrodinger-Newton equation as the defining equa- 
tion for which states could be stable under gravitationally 
induced collapse, and which could not @, 0]- This Schro- 
dinger-Newton equation is a non-linear set of equations 
defined as 
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V 2 U = AirGm 2 |V>| 2 , 



(1) 



where tp is the wavefunction for a particle with mass m, 
H Planck's constant, U the gravitational potential, E the 



energy eigenvalue and G the gravitational constant. Be- 
cause the total energy operator is the generator for time 
translations, we should in principle also be able to use 
these equations to describe the dynamical reduction of 
a quantum superposition. Here we investigate the possi- 
bility of taking the Schrodinger-Newton equation ([l} not 
only as an equation governing the final shape of gravi- 
tationally collapsed states, but also as a generalization 
of the quantum mechanical Hamiltonian operator, which 
can give us both the total energy and the time evolution 
of a quantum superposition state. 

In our study of the possible implied time evolutions, 
we will adopt the philosophy that in order for the 
Schrodinger-Newton equation to be considered as a can- 
didate description for the generator of quantum state re- 
duction, it should at least be able to correctly describe 
the simplest possible situations in which quantum state 
reduction is expected to occur. We will thus calculate 
its implied dynamics in a number of simplified settings, 
without considering whether these situations correspond 
to practical experiments. The rationale is that if the dy- 
namics does not yield the desired outcome in these simple 
cases, we should not expect it to do so in more complex 
real-life situations. 

In the next section we will start out by examinig the 
dynamics of a system defined in a Hilbert space that 
consists of only two possible states. We will see that the 
Schrodinger-Newton equation in that case does not gives 
rise to quantum state reduction, but that it can be eas- 
ily adjusted into a related equation which does give the 
desired result. Next, in section III, we turn to the de- 
scription of a superpostion over three initial states. Here 
even the adjusted equation is seen to harbor only two 
of the three necessary characteristics for a full descrip- 
tion of reduction dynamics. In section IV we will thus 
conclude that although the Schrodinger-Newton equation 
does offer a rationalisation for some of the hitherto un- 
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explained characteristics of quantum state reduction, it 
can nonetheless not be the full story. 



II: A TWO STATE MEASUREMENT 

The very simplest possible experiment that any descrip- 
tion of a quantum state reduction process should be able 
to describe is the evolution of a macroscopic superposi- 
tion over two distinct states to just one of these states. 
One can imagine the macroscopic superposition to be 
formed in a process in which we start out with some mi- 
croscopic superposition of two different quantum states. 
The difference between these states might give rise to a 
gravitational self-energy, but the mass involved is sup- 
posed to be so small that the microscopic matter will not 
collapse by itself. At a certain point in time a coupling 
between the microscopic state and a macroscopic measur- 
ing machine will be instantaneously turned on. In general 
this will yield a macroscopic superposition and the differ- 
ence between the superposed states will generically have 
a finite gravitational self-energy. Now the macroscopic 
mass involved in this self-energy is expected to make the 
collapse process very fast, and as a result it will seem as 
if an instantaneous measurement has yielded only one of 
the two possible outcomes. According to Born's rule Q 
the distribution between the two outcomes found in many 
repetitions of the experiment must mirror the squared 
wavefunction of the microscopic superposition that we 
started out with. 



The General Two State Time Evolution 

To see whether the Schrodinger-Newton equation can 
indeed lead to the collapse of a two state measurement, 
we will first write down the generic time evolution of a 
superposition of two states. The most general superpo- 
sition state over a basis with two elements is given by 
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e^cos(6>/2)|0)+e-^sin((9/2)|l) , (2) 



where n is the norm of the wavefunction, which is usually 
set to 1, and where x is the total phase which is usually 
ignored because it cannot be measured by any quantum 
mechanical process. To define the time evolution of this 
wavefunction we introduce a generator for time transla- 
tions G In its most general form, this generator is 
a complex 2x2 matrix in the basis {|0) , |1)}, with eight 
independent real entries: 



G = 



OLR 
1R 



iai Pr + ipi 
ill Sr + iSi 



(3) 



The Hermitian part of this generator will coincide with 
the usual Hermitian quantum mechanical Hamiltonian 
H. The remaining non-Hcrmitian (and possibly even 



non-linear) part is used to generate the non-unitary col- 
lapse dynamics. The time evolution of the wavefunction 
can be generated from the definition of G by looking at 
the infinitesimal time translation of Vo) 



(1 + leG + O (e 2 )) |Vo> 
e l * cos (6/2) |0) + e -i * sin (6/2) |0)1 . (4) 



Here e is an infinitesimally small parameter that mea- 
sures time in units of H. The parameters TV, X, $ and 
6 defining the new wavefunction at time t — e are writ- 
ten in terms of the old parameters n, Xi f an d at time 
t = 0, and are defined to first order in e. It is now easy to 
extract the time evolution of these parameters using the 
definition of the time derivative. The time derivative of 
for example n is given by the limit e — > of (N — n) / e. 
After some algebra this implies that 

9 = (ai- Si) sin (9) 

+ 2 (fa cos (tp) - (3 R sin (tp)) sin 2 (0/2) 

- 2 (7/ cos (ip) + jr sin (tp)) cos 2 (9/2) 

<P = (ur -Sr) + (Pr cos (tp) + fa sin (tp)) tan (9/2) 

- (jr cos (tp) - 7/ sin (tp)) cot (9/2) 

X = (a R + 5r) + ((3r cos (ip) + Pi sin (ip)) tan (9/2) 

+ (j R cos (tp) - 7/ sin (tp)) cot (9/2) 
11 = - ai cos 2 (6»/2) - (5/ sin 2 (9/2) 
+ 1/2 [fi R - tr) sin (tp) sin (9) 

- 1/2 (p! + ji) cos (tp) sin (9) . (5) 

Notice that in the case of purely unitary time evolution, 
generated by a purely Hermitian generator G, the deriva- 
tives simplify considerably, and become identical to the 
usual quantum mechanical time evolution 



<P 
X 
n 



2 (Pi cos (tp) - Pr sin (tp)) 

Sr) - 2 (Pr cos (tp) + Pi sin (tp)) / tan (9) 
Sr) + 2 (p R cos (tp) + Pi sin (tp)) / sin (9) 
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Specific Time Evolutions 



It is clear from the time derivatives §5§ that the total 
phase and norm variables do not influence the time evolu- 
tion of the superposition state as long as the generator G 
does not explicitly depend on them. We can thus study 
the time evolution by considering only the variables tp 
and 9. Their time evolution can be visualized as a flow 
on the Bloch sphere (see figure Q]). Each flowline on the 
sphere then represents the path traced out by the time 
evolution of an initial state somewhere along the path. 
As long as the time evolution is purely unitary, and the 



FIG. 1: Left: the Bloch sphere. The line m is a line of con- 
stant 9, while I has constant tp. The states |±) are defined 
as \/l/2 [|0) ± Middle: the flow as described by equa- 
tions (J7J- Right: the flow pattern of some general Hamilto- 
nian generator of time evolution. 



FIG. 2: Left: the flow of the type a 7 = -5/ / 0. Right: the 
flow of the type /3r 7^ 0, jr = 0. 



Schrodinger-Newton 



generator G thus purely Hermitian, the flow pattern is in 
fact always the same: it consists of rotations around an 
axis spanned by the two eigenstates of the Hamiltonian, 
as depicted in figure [1] These eigenstates are always on 
opposite poles, and their rotation away from the north 
and south pole depends on the rotation of the Hamilto- 
nian away from being a diagonal operator. 



Hamiltonian Flow 

As an example, let us consider the Hamiltonian with 
/3r = 7.R = A and all other entries zero. 

9 = -2Asin(<p) 

p = — 2A cos (tp) I tan (9) 

X = 2 A cos (if) I sin (9) 

n = 0. (7) 

As mentioned before, we are really only interested in the 
explicit time evolution of tp and 9. To find a closed form 
for the description of the flowlines it is useful to notice 

that dg 



I sin (9) 



-d v [tp sin (9)]. This implies that the 



set of differential equations that we are trying to solve is a 
so called exact set of ordinary differential equations [ll[ , 
and that we can solve it by looking for a potential V 
which obeys 



- dgV = tp sin (9) 
d v V = 0sin(0). 



(8) 



This set of equations is easily solved, and yields the po- 
tential V = 2 A cos(<p) sin(6*). The streamlines describing 
the flow on the Bloch sphere are lines of constant poten- 
tial V, which are given by 



cos (tp) sin (9) — constant. 



(9) 



The flow therefore is a rotation around the axis through 
the north and south poles at \/l/2 [|0) ± |1)], as seen in 
figureU These are of course also precisely the eigenstates 
of the Hamiltonian. The flow in circles around the poles 
(for example the one starting out at |0)) correspods to 
what one usually refers to as Rabbi oscillations. 



In addition to the Hamiltonian part of G we could also 
include non- Hermitian terms in the definition of the gen- 
erator of time evolution. In general these terms can lead 
to many different possible flow patterns on the Bloch 
sphere. To get an idea of the types of patterns that can 
be constructed this way, let us consider three of the most 
basic possibilities for creating purely non-unitary dynam- 
ics. If we set all entries in the 2x2 matrix G equal to zero 
except for the terms ai and Si, the resulting time evo- 
lution does not alter the individual components of the 
superposed state, and only influences the norm of the 
wavefunction. On the other hand, evolutions of the type 
generated by ai = —Sj 7^ and rotations thereof (like 
[3r = — "fR or (3 1 = 71) cause a flow from one pole to the 
opposite pole. That is, one of the poles becomes a source 
and the opposite a sink for the flowlines. Finally, asym- 
metric terms like (3r 7^ 0, = 0, give rise to asymmetric 
flows through a saddle point with flowlines both emerg- 
ing from them and disappearing into them, as shown in 
figure [2 

The specific non-unitary pattern we are interested in 
here, is of course the one given by the non-linear gener- 
ator implied by the Schrodinger-Newton equation ([T]) [^, 
7]. The time evolution generator introduced in the 
Schrodinger-Newton equation for the specific case of a 
superposition of mass over two distinct positions \x = Q) 
and \x = 1) is 



G = 



U(x = 0) A 
A U(x = l) 



(10) 



where the gravitational potential U{x) is given by 
V 2 U(x) = 7 \x)\ . For simplicity we will absorb the 
constant 7 into the norm of the wavefunction, and again 
ignore the dynamics of that norm. We will also take out 
the normal Hamiltonian part of G by setting A = 0. Af- 
ter all, the reduction process should be caused by the 
non-unitary part of the time evolution. 

To solve for the gravitational potential while avoiding 
the infinite self-energy of a point particle, we will consider 
the states labeled by x to represent a mass distribution 
stretched out over an infinite sheet in the y, z-plane, but 
completely localized in the x-direction. The superposi- 
tion thus achieved should be a good description of for 
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|0> |0> 




FIG. 3: Left: the flow as defined by the generator G of equa- 
tion (|10p . with A set to zero. Right: the flow defined by the 
Schrodinger-Newton equation, including the kinetic energy. 



|0> |0> 




n> n> 



FIG. 4: Left: the flow as defined by the generator G' of equa- 
tion (fTB")) . Right: the flow defined by the adjusted generator 
G" of equation (fll)) . 



example a large block of mass which is superposed over a 
distance small compared to its own length. Using appro- 
priate boundary conditions for the gravitational potential 
and completely ignoring the norm of the wavefunction, 
the essential part of the time evolution generator reduces 
to 

Clearly the matrix G is a non- linear operator, because 
it depends on the value of the parameter 9 which de- 
fines the state on which G acts. On the other hand G 
is still akin to a Hermitian matrix in the sense that its 
transpose equals its complex conjugate. The flow pattern 
associated with this generator is easily found to consist 
of circular flowlines around the north and south pole. 
In contrast with the usual Hamiltonian flow though, the 
circulation on the northern hemisphere is in the oppo- 
site direction of its southern counterpart (see figure [3]). 
The gravitational term has thus introduced a division 
between the northern and the southern hemisphere, but 
it has not caused any sinks or sources to appear on the 
Bloch sphere, and it is thus inadequate as a description 
of the dynamical quantum reduction process. Even if 
we reintroduce the mixing parameter A, this will only 
distort the flow lines from their perfectly circular orbits 
and produce some sort of a tennis ball flow pattern as 
depicted in figure [3] However, it does not introduce any 
sources or sinks that would represent the final states in 
a collapse process. 

Alternative Gravitational Terms 

The lack of sources and, more importantly, sinks in 
the flow pattern associated with the Schrodinger-Newton 
equation implies that it cannot be used as a description 
of the dynamical process of quantum state reduction. We 
can however try and introduce sinks into the dynamics by 
slightly altering the Schrodinger-Newton equation. The 
simplest way to do so, as was also suggested by Diosi Q, 



is to turn the gravitational self-energy term into an imag- 
inary potential, so that the equation becomes 

V 2 V> + iUib = Eib 

2m 

V 2 U = AirGm 2 |i/>| 2 . (12) 

The exact physical meaning of the imaginary potential 
term is for the moment unclear, but its mathematical 
form suggests the possibility of referring to it as a dissi- 
pation term induced into quantum mechanics by gravity. 
For now we will not consider the justification of the non- 
unitary dynamics, but merely study the implications of 
having a generator of time evolution given by 

<?=(- ia ?W o y (13) 

The flowlines generated by this matrix all lie along the 
meridian and the flow goes north on the northern hemi- 
sphere, while it goes south on the southern hemisphere 
(see figure 0]). Thus all states starting out above the 
equator will eventually collapse onto the north pole, and 
all states south of the equator find their destination on 
the south pole. If we introduce a normal, real quantum 
mechanical energy into the dynamics as well, then these 
straight flowlines turn into spirals which flow around the 
north-south axis as well as toward one of the poles. 

The dynamics has now two of the three properties ex- 
pected of a working model for quantum state reduction. 
It identifies the states into which a superposition can 
collapse (the poles) by making sure that spatial super- 
positions disappear. At the same time it explains why 
microscopic superpositions can exist while macroscopic 
superpositions are never seen: the spiraling motion for 
microscopic particles is extremely close to perfect circu- 
lar motion because the mass of the particles is small com- 
pared to their usual quantum mechanical potential and 
kinetic energy. On the other hand the gravitational term 
dominates for macroscopic objects, and thus their super- 
position states will be destroyed in a very short time. In 
its present form, however, the time evolution defined by 
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G cannot be used to reproduce the third requirement: 
Born's rule. A state on the northern hemisphere will 
always collapse onto the north pole and never onto the 
south pole. The only way to cure this problem is to in- 
troduce a random variable into the dynamics In a 
description of the state reduction process on the level of 
the Schrodinger-Newton equation the origin of this ran- 
dom variable remains unclear. It could correspond to an 
elusive, as yet unnoticed, field which only manifests it- 
self on the edge between quantum mechanics and gravity. 
An alternative justification could be to identify the ran- 
dom variable with the total phase of the wavefunction. 
We then effectively use the total phase of a quantum me- 
chanical state as a random variable. This seems natural 
as it is an unmeasurable quantity in quantum mechanics. 
After introducing it into the altered Schrodinger-Newton 
equation, the time evolution generator becomes 



G" 



-i[cos(9) + f( X )} 




t [cos (0) + /(*)] 



(14) 



Here / (x) is a function of the random variable \ that will 
be fixed by the requirement that the outcome of the dy- 
namics agree with Born's rule. The flow pattern of this 
adjusted G" is the same as the flow pattern we found 
before for G' , only now the seperatrix between stream- 
ing northward and streaming southward lies not at the 
equator but at the line cos (6>) = —f(x) ( see figured]). 
If we assume x to be taken at random from a flat distri- 
bution between and 2ir, then it is easily checked that 
the dynamics agrees perfectly with Born's rule if we set 

/(x) = xA-i. 

We have thus found a model that (i) describes the col- 
lapse of a quantum mechanical superposition over two 
different states, (ii) distinguishes between microscopic 
and macroscopic superpositions and (Hi) results in the 
emergence of Born's rule if it is repeated many times with 
the same initial conditions, apart from a single random 
variable that remains unobservable. 



Ill: A THREE STATE MEASUREMENT 



Although the model seems to work well for describing 
the reduction dynamics of a two-state measurement, it 
cannot be related to any possible scenario for the solution 
of the measurement problem unless it also works for more 
general superpositions. The first step toward testing the 
model for such a general situation is to ensure that it 
works for a wavefunction superposed over three states 
instead of just two. To do so we repeat the analysis of 





FIG. 5: Left: the quarter of a sphere on which the three- 
state time evolution can be depicted. The line m is a line of 
constant 9 while I has constant rj. Right: the generic flow 
pattern for the flow defined by equations (|17|l . 



the previous section using the initial state 



l^o) 



+ e" 



cos (6»/2) cos (?y/2) |0) 
^ cos (6/2) sin (77/2) |1) 



sin (9/2) \2) 



(15) 



If we also use a 3x3 matrix for the time evolution gener- 
ator, then the computation of the time derivatives of 9, 77 
and so on is exactly analogous to the two state case. To 
see what the effect of our modified Schrodinger-Newton 
time evolution G" is in this case, we consider the states 

0), |1) and 1 2) to represent infinitely high columns of 
mass, positioned on the vertices of an equilateral triangle 
in the x, 77-plane. That way we find that the gravitational 
potential, up to constant prefactor, is given by 

S7 2 U(\x)) cx M x)\ 2 



U(\0)) <x 1 
U(\l)) cx 1 
E/(|2» <x 1 



3 cos 2 (9/2) cos 2 (77/2) 
3 cos 2 (6»/2) sin 2 (r?/2) 
3 sin 2 (9/2) . 



(16) 



Using these values we can then define the generator of 
time evolution G" in analogy with equation (TTJJ by set- 
ting (x\ G" \x) = i (U (\x)) + f x (x))- In analogy to the 
two state time evolution, we have introduced a function / 
which depends on a single random variable and which is 
designed to generate the randomness we need to be able 
to agree with Born's rule. The time derivatives which de- 
scribe the flow of the state vector through configuration 
space during the time evolution defined by this G" are 
then 

9 = sin (9) [a cos 2 (77/2) + b - cos (9) + 

icos 2 (9/2) sin 2 (7,)] 

77 = sin (77) [a - cos 2 (9/2) cos (77)] , (17) 

with a = fo — fi and b = f\ — fa. The relative phases tp 
and </> are constant in time as long as we do not consider 
any off-diagonal elements in G" , and are therefore irrel- 
evant for the reduction process. To visualize the flow we 
can use the surface of one octant of a sphere on which 9 
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measures altitude and r] latitude, so that the states |0), 
|1) and \2) are at the vertices of the surface, as shown in 
figure \5\ The generic flow diagram of the equations (JTTJ) 
has a central source from which all flowlines emanate. 
The flowlines end either at the sinks located at the ver- 
tices of the octant or at a saddle point on one of the edges 
of the surface. From the saddle point the flow continues 
to the vertices again (see figure [5]) . 

Changing the values of a and b corresponds to moving 
the position of the central source over the entire surface 
and at the same time shifting the saddle points along the 
edges. To be precise, the position of the fixed points are 
given in terms of (j),0) coordinates as 



Pa = ( arccos . 

V \a + 2b + 2 

Po-i = (arccos (a) , 0) 

Pi-2 = (tt, arccos (&)) 

P-2-o = (0, arccos (a + b)) 



3a 



, arccos 



1 + 2a + 46 



(18) 



Here Pa is the central source and Pi-j is the saddle point 
on the edge connecting |i) with \j). Clearly the flow 
pattern for the three state superposition fulfills the first 
two requirements for being considered as a description of 
quantum state reduction. The stable points (the sinks) in 
the flow represent precisely the three possible wavefunc- 
tions that do not involve a superposition over gravita- 
tionally distinct states, and that are therefore acceptable 
as possible outcomes of a quantum measurement. The 
time involved in getting to such a stable state is again 
governed by the ratio between gravitational and kinetic 
energy. The microscopic superpositions will thus be able 
to avoid collapse for a very long time, while macroscopic 
superpositions are doomed to collapse within moments 
after their creation. 



Born's Rule 

The challange that is left for the three-state problem is 
to find a function / such that repeated application of the 
measurement model (fT7)> yields Born's rule. As before the 
introduction of some random variable cannot be avoided, 
and again we will try and use just one random variable, 
which could be either an unobserved field or the total 
phase of the wavefunction. 

When proposing an Ansatz for / we should keep in 
mind that the three-state time evolution must reduce to 
the two-state time evolution which we found before in 
the case that the initial state happens to be on one of the 
edges of configuration space. This in fact implies that the 
saddle points on the edges must move along the edges for 
varying \ precisely like the seperatrix moved along the 
meridian on the two-state Bloch sphere. In addition we 
ought to demand that the collection of all possible flow 
patterns posses a 3-fold rotational symmetry in the sense 






FIG. 6: A few of the flow diagrams that are encountered as 
X moves from to 2ir in equations (|19|) . 



that for every flow pattern in the collection there must be 
two more flow patterns which coincide with the original 
one if the vertices are interchanged in a cyclic fashion. 
In the end there seems to be just one possible choice for 
the function / (or equivalently, for a and b) that satisfies 
all of these conditions and depends only on one random 
variable %. This somewhat pathlogical looking choice is 
given by 



-l + §^ 

3x 



1 - 



3X 

-2 + ^ 



if 
if 

if 
if 



X < 

x> 



4tt 
3 

■in 



X> if 



3 
3 



(19) 



These forms for a and b imply that as a function of \ 
the central source moves all around the perimeter of con- 
figuration space, while the saddle points move up and 
down their respective edges, as depicted in figure El This 
ensures that Born's rule will hold on the edges of con- 
figuration space, as it did in the two-state superposition 
scenario. Whether or not it holds away from the edges 
is difficult to prove analytically because an equation for 
the flow lines connecting the central source to the sad- 
dle points is not easily found. Numerically however it 
is rather straightforward to just simulate the reduction 
dynamics many times and compare the result with the 
expected outcome based on Born's rule. 

As it turns out, the proposed dynamics, including the 
definitions (| 19[) . do not agree with Born's rule. The dif- 
ference is shown in figure [7] To fix the mismatch one 
could try other ways to define /. We could look for a 
different scheme in which / satisfies all necessary con- 
ditions but differs from (fT5|) : we could introduce a de- 
pendence of / on i] or 9; or we could introduce addi- 
tional random variables. We did not find any of these 
approaches to be viable. Even the simulated solution in 
which we force the position of the central source to form 
a flat distribution over configuration space in the course 
of many experiments does not yield the desired result 
(see figure [7]). Moreover, the introduction of more ran- 
dom variables would be a rather undesirable element in 
the theory, because they should physically emerge from 
ever more fluctuating fields for whose existence there is 
no empirical evidence. 
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FIG. 7: False color plot of the difference between the probability for finding outcome |0), following the time evolution G" from 
the initial state (r/,6) and the corresponding probability according to Born's rule (i.e. cos 2 (#/2) cos 2 (?7/2)). 
Left: Using the measurement scheme defined by equations (|19[) . in which the source of the flow wanders around all edges. 
Right: Using the measurement scheme in which we force the centre of the flow to be in a random position for each value of \- 



IV: CONCLUSIONS 



The results discussed in this paper show that the 
Schrodinger-Newton equation ([T]) which was proposed by 
Penrose as a replacement for the quantum mechanical 
Hamiltonian @, 0] can hardly serve as a description of 
the dynamics of quantum state reduction. On the other 
hand, a slight modification of the equation, i.e. making 
the gravitational potential energy appear as an imaginary 
term, causes the associated dynamics to show at least two 
of the three characteristics necessary for being considered 
as a possible reduction model. The modified equation 
causes the system to evolve toward states which are not 
superposed over gravitationally different positions, and 
thus selects the correct Pointer basis for the quantum sys- 
tem to collapse into(§. H^|. On top of this the equation 
naturally provides a reason for the observed difference 
between microscopic and macroscopic objects. Micro- 
scopic systems have a very small gravitational potential 
energy as compared to their internal quantum mechani- 
cal potential and kinetic energy. The collapse process will 
therefore be so slow that it cannot be noticed on human 
timescales. On the other hand the gravitational term will 
dominate in the dynamics of macroscopic superposition 
states, and these will thus collapse before their existence 
can be noticed 0, Q • 

If the wavefunction is a sum of only two distinct states 
then the addition of a random variable into the dynamics 
rather straightforwardly leads to the desired statistics for 
the outcomes of measurements. However, as soon as the 
wavefunction represents a superposition over more then 
two states, it becomes impossible to force the dynam- 
ics of the collapse model to agree with Born's rule using 
only one random variable. Even apart from the fact that 



there is no physical ground for introducing them, more 
than one random variable does not automatically solve 
the problem. At least the obvious choices of how to im- 
plement them into the theory do not seem to yield the 
desired outcome. 

We therefore conclude that the Schrodinger-Newton 
equation, even in its modified forms, is still not fit as 
a complete description of the dynamical collapse process 
of quantum superpositions, because it is unclear how to 
make it agree with Born's rule under any but the most 
basic circumstances. 
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